Non-Markovian decoherence of localized nanotube excitons by acoustic phonons 
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We demonstrate that electron-phonon interaction in quantum dots embedded in one-dimensional systems 
leads to pronounced, non-Markovian decoherence of optical transitions. The experiments we present focus on 
the line shape of photoluminescence from low-temperature axially localized carbon nanotube excitons. The 
independent boson model that we use to model the phonon interactions reproduces with very high accuracy the 
broad and asymmetric emission lines and the weak red-detuned radial breathing mode replicas observed in the 
experiments. The intrinsic phonon-induced pure dephasing of the zero-phonon line is two orders of magnitude 
larger than the lifetime broadening and is a hallmark of the reduced dimensionality of the phonon bath. The 
non-Markovian nature of this decoherence mechanism may have adverse consequences for applications of one- 
dimensional systems in quantum information processing. 



In single-wall carbon nanotubes (SWNTs), the one- 
dimensional (ID) nature of electronic states results in an en- 
hancement of Coulomb interactions and has lead to the ob- 
servation of remarkable effects such as Luttinger-liquid states 
in transport experiments [2]. On the other hand, a common 
feature of the small diameter semiconducting SWNTs is the 
unintentional confinement of the electronic or excitonic states 
along the nanotube axis; this leads to the formation of zero- 
dimensional (OD) quantum dot (QD)-like states, manifesting 
themselves through Coulomb blockade in transport experi- 
ments A 01 or through localized emission |5] and strong 
anti-bunching in photoluminescence (PL) [6]. Spatial confine- 
ment of electrons or excitons in SWNTs is in fact desirable 
for applications in single photonics or quantum information 
processing (QIP), which in turn motivates theoretical and ex- 
perimental investigations of elementary physical properties of 
optically active QDs in SWNTs. 

Electron coupling to Raman-active phonon modes such as 
radial breathing mode (RBM) and G mode has been an in- 
valuable tool for characterizing SWNTs. Interest in long- 
wavelength acoustic phonons of SWNTs primarily stems from 
their potential use as nano-mechanical resonators. A question 
that is of interest from QIP perspective is the decoherence in- 
duced by interactions between QD electrons or excitons and 
ID acoustic phonons: while this problem has been addressed 
theoretically for QD excitons in quasi- ID III-V nanowires in 
the high temperature limit 0] experimental confirmation of 
predicted features have not been obtained. 

In this letter, we analyze the influence of exciton-phonon in- 
teractions on the PL line shape of SWNT QDs. We show that 
the experimental PL spectra consisting of broad and strongly 
asymmetric lines can be fit with remarkable accuracy using 
an analytical solution to the independent boson model, over 
a dynamical range spanning two orders of magnitude and a 
temperature range varying from 5K to 35K. In our model, we 
assume that the lowest bright SWNT excitons, formed out of 
bound electron-hole pairs from the first band-to-band transi- 
tion En |H1l are confined to form QD-like states; this assump- 
tion is justified since PL has been shown to exhibit strong pho- 
ton antibunching [6]. The principal result of our work is the 
demonstration that ID nature of acoustic phonons and the re- 



sulting Ohmic exciton-phonon coupling gives rise to strong 
pure dephasing of spatially localized optical excitations even 
as we approach the zero-temperature limit. Even though we 
study localized SWNT excitons in our experiments, we em- 
phasize that our conclusions are valid for localized transitions 
in any ID system. Our findings also suggest that simple PL 
measurements can be used to determine quantities such as the 
localization length scale of excitons and the deformation po- 
tential coupling strength in carbon nanotubes. 

We measured the PL from isolated SWNTs at cryogenic 
temperatures. Details on the experimental setup and sample 
preparation are given in jfj. Briefly, we deposed surfactant- 
embedded CoMoCat SWNTs directly on a functionalized 
solid-immersion lens (SIL) and used a confocal microscope 
to detect PL upon excitation with a tunable TkSapphire laser 
(CW or pulsed mode). The sample is placed in a helium bath 
cryostat and the temperature is controlled locally through a 
thermo-resistor. As discussed in [6], spatially localized emis- 
sion and strong photon anti-bunching indicate that the PL 
signal originates from unintentionally axially confined exci- 
tons, i.e. from SWNT QDs. Figure 1 shows typical experi- 
mental PL spectra (open circles) at T = 4.2 K that we ob- 
served: the line shape is clearly asymmetric with a full-width- 
half-maximum (FWHM) of 3.5 meV. Time-resolved measure- 
ments indicate that the PL decay follows a bi-exponential 
law with a fastest time constant of 36 psec (Fig. 1 inset) in 
agreement with previous reports (see JH [Io|] and references 
therein), which implies that the measured linewidth is a fac- 
tor of 100 larger than the lifetime broadening. This partic- 
ular nanotube PL exhibited strong photon antibunching with 
a value for the normalized second-order correlation function 

s (2) (o)=o.ia. 

To explain this broad and asymmetric PL line shape we de- 
velop a theoretical model where we assume a harmonic con- 
finement potential for excitons along the nanotube axis. For 
confinement length scales that are long compared to the ex- 
citon Bohr radius, QD states are superpositions of free exci- 
tonic states with envelope wave functions determined by the 
ID harmonic oscillator Hamiltonian. We further assume that 
the confinement is strong enough to lead to well-separated QD 
levels, with an energy splitting AE ^S> ksT, where fc^ is the 
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FIG. 1: Photoluminescence of a nanotube quantum dot at T = 4.2 
K and P = 300 nW (open circles) and the calculated spectrum (solid 
line) plotted on a logarithmic scale. For the theoretical model we 
used the parameters Ds = 14 eV, Drbm = 1-4 eV/A, Tiujrbm = 
37.2 meV, AE = 44 meV and T fit = 6.3 K. The upper and lower 
blue solid lines are calculated spectra with identical parameters but 
Drbm = 1-8 eV/A and Drbm = 1-0 eV/A, respectively, demon- 
strating the precision with which the coupling strength can be de- 
termined. The inset shows a time-resolved measurement (dark line) 
The red line is a bi-exponential fit of the PL decay revealing a main 
fast component with characteristic time t\ = 36 ± 0.2 ps and weight 
Ai = 100 ± 0.3 and a vanishingly weak long tail (t 2 = 284± 40 ps 
, Ai = 1.2 ± 0.2). From the value of t\ we expect a lifetime-limited 
line-width of ~ 0.036 meV, in strong contrast to the measured and 
calculated width of 3.5 meV. 



Boltzmann constant and T the absolute temperature. 

To account for the coupling of this QD to the acoustic 
phonons of the SWNT, we follow the extension of the inde- 
pendent boson model developed in [11]. We assume that the 
phonon dispersion of the bare SWNT remains unaltered by 
the confining potential. Since no atomic bonds are formed 
with the embedding surfactant, we also assume that coupling 
to phonons of the surrounding material is negligible. Long 
wavelength acoustic phonons in nanotubes are well described 
by a continuum elastic model J72 , 13 ] and they couple to elec- 

Il4ll . Because of 



trons through deformation potentials [12, 13, 
momentum conservation, only phonon modes with zero an- 
gular momentum can couple to excitonic states within a sin- 
gle subband En; the potentially relevant modes are therefore 
stretching, twisting and radial breathing modes. We neglect 
the coupling to the twisting mode, whose strength is expected 
to be one order of magnitude smaller than to the stretching 
mode [12]. The stretching mode corresponds to lattice vibra- 
tions with longitudinal displacements of the atoms along the 
nanotube axis yjfl . It has linear dispersion ll> s (q) = v s q for 
small wave- vectors lfl3ll and the strength of its coupling to ex- 
citons is characterized by the deformation potential Ds- It 
would provide the dominant contribution to the broadening 
of the localized exciton emission, a.k.a. the zero-phonon line 
(ZPL). Since the RBM has an approximately flat dispersion 
for small wave-vectors (with an energy Hujubm), its principal 
contribution would be the appearance of phonon sidebands in 



PL that are separated from the exciton line by friORBM- 

Following ref 111 111 , we obtain the exciton-phonon matrix 
elements: gj(q) — Gj(q) ■ T exc {q) where j stands for S 
(stretching mode) or RBM. The form factor !F exc (q) = 
J dz\ i ^ exc (z)\ 2 e lq ' z provides a cut-off for the coupling in mo- 
mentum space that is inversely proportional to the spatial ex- 
tension of the wave function. With our choice of a parabolic 
potential, we get for the QD ground state a gaussian envelope 
along z: fy exc (z) = ^1/4^1/2 ex P~^r where a is the con- 

c (q) = exp--^ and the 
where m* _ is the effective 
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finement length. This yields T' 
energy splitting: AE = 



mass of the exciton. The deformation potential couplings are: 

Q s (q) = . D '1 and Q R BM{q) = ,®A BM where L 

is the length of the nanotube and p its linear mass density. 

The linear susceptibility x(t) (t > 0) of the QD in response 
to a f5-shaped laser pulse at t = |n , 15 ] can be decomposed 
into a temperature-dependent and a temperature-independent 
contribution = —ie~ int XT{t) ■ Xo(t) with 

XT (t) ex iexpJ2hs(l)\ 2 [-ns(q)\e-^ q)t -l' 
1 
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where jj (q) 

( e hu 3 (q)/k B T _ -Q 

number. The 

n = n - 



9j{.l)/ UJ ]{q)^ and n ](l) = 
is the phonon occupation 
polaron-shifted transition energy is 
\jj(q)\ 2 ujj(q) , where Ml is the bare 
energy of the QD-state. Performing the Fourier transform 
of x(t) an d taking its imaginary part gives the absorption 
spectrum of the QD. The PL-emission line shape that we 
measure experimentally is obtained by taking the mirror 
image of the absorption profile with respect to the ZPL II 1611 . 

The solid curve in Fig. 1 shows that this simple theoreti- 
cal model describes the experimental observation remarkably 
well. To obtain this fit, we use SWNT parameters that were 
reported earlier; in particular, we take p = 1.67 • 10~ 15 kg/m, 



v s = 19.9 Km/s 1 12] and to* 



0.2 



[17]. An electron- 



phonon deformation potential for the stretching mode [ 12, 13] 
of strength Ds = 14 eV has been extracted from data on low- 
field mobilities in semiconducting SWNTs |14] j26ll . Our pro- 
cedure, applied to 10 different SWNT spectra, involves vary- 
ing the values of Ds and the confinement to obtain the best 
fit to the data. The range of Ds values we thereby obtained 
varied from 12 to 14 eV. The other fitting parameter is the con- 
finement potential along the SWNT axis which determines the 
cut-off in momentum space, or equivalently the level spacing 
AE. We found values for AE between 20 and 44 meV, con- 
sistent with our assumption AE ^> ksT w 2.8 meV (at 35 K) 
and corresponding to an extension of the wave function enve- 
lope on the order of ~ 10 nm. Since this confinement length- 
scale is bigger than the Bohr radius [8|,ll8|]> our assumption of 
exciton center-of-mass confinement is also justified. 

The PL spectrum from Fig. 1 is plotted on a logarithmic ver- 
tical scale to show the quality of the fit over two decades of 
counts and to emphasize the presence of a red-detuned replica 
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FIG. 2: Left panel: Experimental (open circles) and calculated (red 
solid lines) PL spectra for the same SWNT at different temperatures 
T mes (the vertical scale is arbitrary). The excitation power was kept 
constant at P — 3fiW. The effective temperature used in the sim- 
ulation is given by Tfu, the other parameters are Ds = 12.7 eV 
and AE — 22 meV. Right panel: Emission from another nanotube 
at constant sample temperature T = 4.2 K. The changes induced by 
increasing the laser intensity from 4 nW to 40 /iW are very well 
reproduced by increasing the temperature Tju in the model. The 
parameters for this SWNT are D s = 13 eV and AE = 20 meV . 



that is exactly reproduced by the model. The RBM deforma- 
tion potential for excitons are well documented 11191 |20L 12111 
and expected to lie between 1 and 2 eV/A. Here we obtained 
the best fit for D RBM = 1.4 ± 0.2 eV/A and a RBM en- 
ergy of 37.2 meV corresponding to a SWNT diameter of about 
0.7 nm, in good agreement with CoMoCat SWNTs emitting 
at this wavelength. The linear power dependence of the red- 
detuned replica for very low excitation intensities (not shown) 
is consistent with that of a phonon sideband and excludes an 
explanation based on biexciton emission. 

Figure 2 shows the experimental PL lines (open circles) for 
two different SWNTs. In the left panel, the temperature is in- 
creased by flowing current in the thermo-resistor. The very 
asymmetric lines at the lowest temperatures and the evolution 
to broader symmetric lines with increasing temperature is a 
general feature common to virtually all of the 30+ SWNTs 
that we measured. The computed spectra (solid red lines) 
reproduce both qualitatively and quantitatively this behavior 
when the coupling to the stretching mode only is taken into 
account. Only the considerable shifts in the transition ener- 
gies are not reproduced by the model; arbitrary offsets were 
therefore introduced for each spectrum. This unexpected sen- 
sitivity on minor environmental changes is another evidence 
for localized emitting states [fj. In contrast, the transition en- 
ergy of a delocalized exciton is not expected to shift so dra- 
matically with temperature j22Tl . 



FIG. 3: Upper panel: Calculated emission line-shapes at T — 5 K 
for a QD coupled to a phonon bath of dimension d — 1 (left) and 
d — 3 (right) are shown. The coupling to a single acoustic phonon 
branch with linear dispersion (v s = 19.9 Km/s) is considered. To 
extend the model to higher dimension, we assume an isotropic cou- 
pling of the QD to LA phonons described by the same deformation 
potential. A lifetime broadening corresponding to ti = 36 ps is in- 
troduced. Lower panel: For illustration, fit of the spectrum shown 
in Fig. 1 (open circles) with only the temperature dependent (Eq. Q] 
red line) and the temperature-independent (Eq. [2] blue line) contri- 
butions. Each curve is scaled to span the same total area. The inset 
shows the susceptibility (T = 6 K) used to fit the data (black curve), 
where a deviation from an exponential decay is visible for t < lpsec. 
The blue curve of the inset shows the same function at T = 0.1 K 
where the non-Markovian nature of dephasing becomes clearer. 



The right panel in Fig. 2 (different SWNT) shows that the 
effect of increasing the excitation power can be reproduced 
by only increasing the temperature of the phonon bath in the 
model and keeping all the other parameters fixed. From this, 
one can conclude that the laser excitation led a to a very local 
heating of the SWNT that was not directly measured by our 
thermo-resistor placed about 1 mm away from the nanotubes. 

We emphasize that the unusually broad lines observed in 
SWNT PL can be accurately and quantitatively explained by 
solely taking into account interaction with long wavelength 
ID acoustic phonons. The resulting ohmic coupling [23] gives 
rise to linewidths that are two orders of magnitude larger than 
the ones expected from lifetime broadening. To highlight this 
feature, we have plotted in Fig. 3 a comparison of the theo- 
retical PL line shape for a QD interacting with ID and 3D 
phonon reservoirs. Coupling to 2D or 3D acoustic phonons 
yields a superohmic spectral function 12311 . for which the in- 
dependent boson model predicts a (lifetime-broadened) ZPL 
accompanied by phonon sidebands (Fig. 3, upper right). This 
is in strong contrast to our observations from SWNTs. M27I1 

The enhanced pure dephasing of QD transitions in ID struc- 
tures was analyzed theoretically in Ref. Q7[| where the fi- 
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nite width was attributed to the high temperature (Markovian) 
limit. This part is due to the temperature-dependent contri- 
bution xt(0 m Eq.(Q~|i and vanishes linearly as T — > 0. It 
provides a symmetric contribution to the line shape (as illus- 
trated in Fig. [3] lower panel, red curve). In contrast, the highly 
asymmetric character of the observed line shapes at low tem- 
peratures is largely due to the temperature-independent part of 
the susceptibility, Xo(t), in Eq.©. For comparison, we plot 
these two contributions separately, although they appear in the 
line shape as a convolution. The temperature-dependent term 
is responsible for the presence of a high-energy tail at non- 
zero temperature, while the broad low-energy tail we measure 
constitutes a direct signature of the temperature-independent 
pure dephasing. Furthermore, we find that the characteristic 
energy of the latter quasi-exponential tail is directly propor- 
tional to the width of the form factor in g-space. This depen- 
dence is to be expected since the inverse localization length 
scale determines the highest phonon energies to which the QD 
is effectively coupled. A stronger confinement therefore leads 
to a faster dephasing, even at arbitrarily low temperatures. 

In Fig. 3 (lower panel) the inset shows the time dependence 
of the calculated linear susceptibility x(t) use d to fit the ex- 
perimental spectrum: the deviation from an exponential decay 
that is apparent for t < 1 psec is a hallmark of non-Markovian 
decoherence B24I1 . Our calculations indicate that this initial 
(Gaussian) decay crosses over to a power-law decay and fi- 
nally to an exponential decay at long times. The regime where 
power-law decay is dominant expands as temperature is de- 
creased; further reduction of the temperature to 0. IK leads to 
the appearance of power-law decay of the susceptibility ex- 
tending to 15 psec (blue curve). Observation of such a strong 
deviation from Markovian dynamics should be within reach. 

To conclude, we have used an extension of the indepen- 
dent boson model to reproduce with remarkable quantitative 
accuracy the asymmetric and broad lines that we measured 
on isolated SWNTs. The model also reveals the changes in 
the shape and width with temperature. In strong contrast to 
3D systems, the reduced dimensionality of the phonon bath 
in SWNTs leads to ohmic coupling of confined excitons to 
stretching mode phonons: in this limit, the ZPL is replaced 
by a power-law singularity. This conclusion remains valid at 
arbitrarily low temperatures, where the induced pure dephas- 
ing becomes manifestly non-Markovian, and would apply to 
QDs defined in any one-dimensional system. These intrinsic 



effects are of fundamental interest in view of applications in 
quantum information processing, where non-Markovian na- 
ture of decoherence may substantially alter the attainable ac- 
curacy thresholds M25I1 . 
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